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Abstract 

The Zamolodchikov model describes an exact relativistic factorized scattering theory of 
straight strings in (2 + l)-dimcnsional space-time. It also defines an integrable 3D lattice 
model of statistical mechanics and quantum field theory. The three-string S-matrix satisfies 
the tetrahedron equation which is a 3D analog of the Yang-Baxter equation. Each S'-matrix 
depends on three dihedral angles formed by three intersecting planes, whereas the tetra- 
hedron equation contains five independent spectral parameters, associated with angles of 
an Euclidean tetrahedron. The vertex weights are given by rather complicated expressions 
involving square roots of trigonometric function of the spectral parameters, which is quite 
unusual from the point of view of 2D solvable lattice models. In this paper we consider a 
particular four-parameter specialization of the tetrahedron equation when one of its vertices 
goes to infinity and the tetrahedron itself degenerates into an infinite prism. We show that 
in this limit all the vertex weights in the tetrahedron equation can be represented as mero- 
morphic functions on an elliptic curve. Moreover we show that a special reduction of the 
tetrahedron equation in this case leads precisely to an example of the tetrahedral Zamolod- 
chikov algebra, previously constructed by Korepanov. This algebra plays important role 
for a "layered" construction of the Shastry's i?-matrix and the 2D S'-matrix appearing in 
the problem of the ADS/CFT correspondence for J\f = 4 SUSY Yang-Mills theory in four 
dimensions. Possible applications of our results in this field are briefly discussed. 



1 Introduction 



The tetrahedron equation [1,2] is a three-dimensional analog of the Yang-Baxter equation. It 
implies the commutativity of layer-to-layer transfer matrices [3] for three-dimensional lattice 
models of statistical mechanics and field theory and, thus, generalizes the most fundamental 
integrability structure of exactly solvable models in two dimensions [4] . 

The first solution of the tetrahedron equation was proposed by Zamolodchikov [1,2] in the 
context of an exact relativistic factorized scattering theory of "straight strings" in 2+1 di- 
mensional space-time. Subsequently, Baxter [5] proved that this solution indeed satisfies the 
tetrahedron equations and then [6] exactly calculated the free energy of the corresponding solv- 
able three-dimensional model in the limit of an infinite lattice. Some further developments in 
the field of 3D integrability in a more general setting could be learned from [7-12]. Many special 
properties of the Zamolodchikov model were studied in [13-16]. In this paper we will consider 
some additional analytic and algebraic problems related to this model. 

The spin variables in the Zamolodchikov model take two values. In the original formulation 
[1, 2] the spins are assigned to faces of the lattice. Baxter [6] used a different, but closely 
related "interaction-round-a-cube" (IRC) approach, where the spins are assigned to 3D cells of 
the lattice, or, equivalently, vertices of the dual lattice. Here we will use yet another approach, 
namely the vertex formulation with edge spins, found by Sergeev, Mangazeev and Stroganov [17], 
which is related to the other formulations by the "cube- vertex" correspondence [18]. The most 
general edge-spin solution of the tetrahedron equation in the model contains five independent 
spectral parameters associated with angles of an Euclidean tetrahedron. Two different reductions 
of this solution with a fewer number of parameters, known as the "static" and "planar" limits, 
were previously found by Korepanov [19] and Hietarinta [20]. 

In this paper we consider another particular limit when one of the tetrahedron vertices 
goes to infinity, whereas the tetrahedron turns into an infinite prism. We show that the all 
Boltzmann weights in this case can be parametrised in terms of meromorphic functions on an 
elliptic curve. Moreover we show that a special reduction of the tetrahedron equation in this 
case leads precisely to the example of the Zamolodchikov tetrahedral algebra, constructed by 
Korepanov [19]. This algebra plays an important role for a "layered" construction [21-23] of the 
Shastry's i?-matrix [24] underlying the integrability of the ID Hubbard model. Remarkably, the 
same i?-matrix [25-28] appears in the problem of the AdS/CFT correspondence for N = 4 SUSY 
Yang-Mills theory in four dimensions. In conclusion we briefly discuss possible applications of 
our results to this field. 

The organization of the paper is as follows. In Sect. 2 we briefly review known results on the 
tetrahedron equation in the Zamolodchikov model, used in this paper. In Sect. 3 we derive an 
elliptic parameterisation of the Boltzmann weights. In Sect. 4 we consider special limits of the 
tetrahedron equation and the tetrahedral Zamolodchikov algebra. 

2 The vertex formulation of the Zamolodchikov model 

The Zamolodchikov model can be formulated on any 3D lattice formed by an arbitrary set of 
intersecting 2D planes in the Euclidean space M 3 , such that no four planes intersect at one 
point. In particular, this definition includes a regular cubic lattice as the simplest possibility. 
The fluctuating spin variables in the model take two values, denoted below and 1. Here 
we will use the vertex formulation of the model found in [17]. The spins in this case are 
assigned to edges of the lattice while local Boltzmann weights are assigned to vertices at three- 
plane intersection points. Let -R^f 2 2 '/ 3 3 denote the weight corresponding to a configuration of six 
edge spins ii,«2)*3,Ji>j2)j3 = 0,1, arranged as in Fig.l. This quantity could be conveniently 
associated with a linear operator 

R : C 2 ®C 2 ®C 2 ^C 2 ®C 2 ®C 2 , (1) 
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Figure 1: Arrangement of spins around an elementary vertex formed by three intersecting planes with 
the dihedral angles 81,62, 83 and the edge angles 01, a<i, as. 



acting in a direct product of three two-dimensional vector spaces, where the pairs of indices 
Ji)> (*2j J2) and (23, J3) serve as matrix indices in the first, second and third spaces, respec- 
tively. The edge states "0" and "1" correspond to the basis vectors vo and v\, 

u ° =(("))' Vl= (l)' u o,«iGC 2 . (2) 

The vertex weights depends on three spectral parameters, 61,62,63, which are identified with 
dihedral angles between the three planes forming the vertex, as shown in Fig.l. To indicate 
this dependence we will write the weights as R(#i, #2, #3)- Consider a spherical triangle with 
the angles 61,62,63 and let 01,02,03 denote three sides of this triangle opposite to the angles 
61,62,63, which obey the spherical sine theorem 

sin6>i sin 6*2 sin6> 3 

A = = = . (3) 

sin ai sin 02 sin 03 

Define related variables 

2 a = 61+62 + 63-71, ai = 6i - a , 

2ft = 2ir - ai - a 2 - a 3 , ft = it - m - ft, 

for i = 1,2,3. They satisfy the relations 

a + a± + a 2 + a 3 = tt, ft + ft + ft + ft = ^ • (5) 

The variables ao, a.\, a 2 , 013 are called excesses of the spherical triangle. Choose 61,62, 63 so that 
61,62, 63, ao, . . . ,03 and ft, . . . , ft are all real and between and n. Further, define four 
variables (taking positive values of the square roots) 



(4) 



\ tan 'J. i = 0,l,2,3, (6) 



3 



which are constrained by the relation 

1 — tytf — t\t\ ~~ *0*3 ~~ *1*2 ~~ ^1*3 ~~ *2*3 + ^0*1*2*1 = 0' (7) 

following from (5). The weights obey a "parity conservation law", such that 32 vertex configu- 
rations with an odd number of 0's have the identically vanishing weight 

Kith = °> if + ji + h + 33 (mod 2) . (8) 

The remaining 32 non-vanishing matrix elements of R(#i, 02,0s) have the form [17] 1 



R 0,0,0 
-"0,0,0 


_ R 0,1,1 _ 

- -"0,1,1 - 


pl,0,l _ 

-"1,0,1 - 


pl,l,0 

-"1,1,0 


1 

— 1 


111 

-"1,1,1 


1 n n 
_ pi, 0,0 _ 

— -"-1,0,0 — 


n i n 
p0,l,0 _ 

-"0,1,0 — 


n n 1 
p0,0,l 

-"0,0,1 


= t()tit 2 ts 


o0,l,0 
-"-0,0,1 


R 0,0,l 
— -"-0,1,0 — 


pl,0,0 _ 

-^1,1,1 - 


pi, 1,1 

-"1,0,0 


= t 2 ts 


ol,0,l 
1,1,0 


pi, 1,0 _ 
~~ -"-1,0,1 ~~ 


p0,l,l _ 

-"0,0,0 — 


p0,0,0 
0,1,1 


= toti 


7,1,1,1 

-"-0,1,0 


_ pl,0,0 _ 
— -"-0,0,1 — 


p0,0,l _ 

-"1,0,0 — 


p0,l,0 

-"1,1,1 


= — i*i*3 


R 0,0,0 
1,0,1 


_ p0,l,l _ 
— -"-1,1,0 ~~ 


pl,l,0 _ 

-"0,1,1 - 


pl,0,l 

-"0,0,0 


= i*o*2 


p0,0,l 

"1,1,1 


— -"-0,0,1 — 


pl,0,0 _ 

-"0,1,0 — 


p0,l,0 

-"1,0,0 


= ht 2 


ol,l,0 
-"-0,0,0 


_ p0,0,0 _ 

- -"1,1,0 - 


p0,l,l _ 

-"1,0,1 - 


pl,0,l 

-"0,1,1 


= *0*3 



(10) 



In the following we will regard the weights as functions of the three independent variables 
0i,02,03, or, equivalently, of the four dependent variables ao,ai,a 2 ,as, constrained by the 
relation (5). We would like to stress that the expressions (10) do not contains any other free 
parameters. 

In the original face-spin formulation of ref. [2] the weights are explicitly symmetric with 
respect to all spatial symmetry transformations, generated by permutations of the lines 1, 2, 
3 and inversions of their directions. Namely, the weights remains unchanged if the associated 
permutations of spins are accompanied by the corresponding transformations of the dihedral 
angles. Contrary, to this the operator H(9i,0 2 ,0s) is not explicitly symmetric. Nevertheless 
it possesses the same spatial symmetry group of the 3D cube [17,29], but its realisation now 
involves linear similarity transformations of the weights. This 48-element group is generated 
just by two transformations, which we choose as 

R(7r-6>i,7r-6»2,0 3 ) = (a y ® a y ® a z ) R tlt2 (0 1 ,6 2 ,6 3 ) (a y ® a y ® a z ), 

(11) 

Pis R* 1 * 2 * 3 (6>3,6> 2 ,#i) Pis = {T®T®T)R{e u 6 2 ,0s) (t®t®t)- 1 , 

where the superscripts ij, i = 1,2,3, denote the matrix transposition in the i-th space, P13 
denotes the operator permuting the spaces 1 and 3 and a x , a y , a z denote the Pauli matrices, 



1 In [17] the weights are given in terms of the variables /3o , /3i , P2 , fi'i ■ Here they are re-expressed in terms of 
a's with the help of the identity 

ai atj f3k , Pi 

tan — tan — tan — tan — (9) 

where (i, j, k, I) is any permutation of the indices (0, 1, 2, 3). 
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In the following the symmetry properties will not play any essential role, since the elliptic 
parameterisation considered below explicitly breaks the spatial symmetry. 
The operator R satisfies the tetrahedron equation, 



R-123(#l, #2, #3) R-145(01> #4) #5) R-246(y ~ #2, #4j #6) R-356(#3> 7T — #5, ^) 



R 



■356 ^3> 



T 1 " — 0&) R246( 7 '" — #2, #4) #6) R-145(#l> #4> #5) R-123 (#1 1 $2 , #3) > 



(13) 



which ensures the integrability of the model. The above equation involves operators acting in 
a direct product V = (C 2 )® 6 , of six identical two-dimensional vector spaces 2 . Geometrically, 
it is associated with an arbitrary intersection of four planes. The indices in (13) numerate six 
intersection lines among these planes, which form extended edges of an Euclidean tetrahedron. 
Each of the operators Ryfc depends on its own set of spectral parameters, determined by the 
dihedral angles corresponding to the lines i, j and k. The geometric arrangement of these angles 
is shown in Fig.2 (note that it is exactly the same as in [17], but different from that used in 
Eq. (2.2) in [5]). Altogether there are six dihedral angles corresponding to six edges of the 





Figure 2: Graphical representation of the tetrahedron equation (13) showing the arrangement of the 

internal dihedral angles of the tetrahedron 



tetrahedron. At this point it is worth mentioning that these angles are not independent; they 
satisfy one nondinear constraint, 



det 



/ 1 

cos 6*i 

COS 02 
\ —COS #4 



cos 6*i 

1 

COS #3 
COS #5 



COS 6*2 
COS 6*3 

1 

cos 9q 



- COS 6*4 \ 
COS #5 

cos 9q 

1 / 



= 0, 



(14) 



which follows from the vanishing of the Gram determinant of four unit normal vectors to faces 
of a tetrahedron in the Euclidean 3-space. Therefore, Eq.(13) contains only five independent 



We follow usual notations, where R123 acts non-trivially in the first three spaces and coincides with the identity 
operator in all the remaining components of the product. The other operators in (13) are defined similarly. 
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parameters. In matrix form it reads 

£ R(9 1 ,9 2 , R ^ - **, *4, ^6)£tt ^3, vr - * 8 )222 

fc 4 fc 5 feg 

= ^ Ufa, VT - 6 )££* ^(VT - 6 2 , 4 , *(«M4, 05)*tf 5 fc, • 

(15) 

Altogether there are 2 12 distinct equations, half of which is non-trivial (for the other half both 
sides vanish identically due to the parity conservation (8)). There are several proofs that the 
weights (10) indeed satisfy these equations. Firstly, thanks to the "cube-vertex" correspondence 
of ref. [18], this follows from the Baxter's proof [5] of the tetrahedron equation for the interaction- 
round-a-cube (IRC) formulation of the model. Baxter used various symmetry relations and 
a computer-aided analysis of a pattern of signs in some equations to reduce all non-trivial 
equations to just two equations, which were then proven with the help of spherical trigonometry. 
Secondly, there exists a completely algebraic proof of (15) given in [17]. This proof (as well as 
its earlier variant [9] for the IRC formulation) is based on a repeated application of the so-called 
restricted star-triangle relation, invented in [8]. It is worth mentioning, that this relation was 
later interpreted as the five-term identity for a cyclic version [30] of the quantum dilogarithm [31]. 
Finally, an alternative (and, in fact, simpler) algebraic proof of (15), based on the auxiliary linear 
problem, was obtained in [32]. 

To conclude this Section, mention a simple fact that Eq.(13) is unaffected by arbitrary linear 
similarity transformations in any of the six vector spaces therein, 



R 



ijk 



G fc ) R ijk (G< <8> Gj <g> G fc ) 



(16) 



where G, 



1, 2, . . . , 6 are arbitrary non-degenerate matrices. We will use this freedom to 



bring the expressions for matrix elements of R(0i,02,03) to the most convenient form. 



3 Elliptic parameter isat ion of the weights 

3.1 Linear transformations of the weights 

For the following analysis it is convenient to introduce the operator 

L123 = (D(£) ® F ® F) R123 (D(0®F®F)-\ 

where 



(17) 
(18) 



which differs from R123 by a simple similarity transformation of the type (16). The parameter £ 
will be specified later on. Evidently, the transformation breaks the symmetry between the three 
spaces in (17), since the first space is distinguished from the other two. The matrix elements of 
L123 can be arranged into a set of 4 x 4 matrices Lq, Lq, L\ and L\, acting in the tensor product 
of the second and third spaces, while the indices 0, 1, labeling these matrices, refer to the first 
space, 



(■^0)2,3 



where 



/ a al \ 

b c 

c b 

\ d a J 



1 — t$t\ + tits + 40*1*2*3) 
1 + *0*1 + *2*3 — *0*1*2*3 



(A: 



2,3 

b 
d 



I b 






1 + *o*i 
1 - *0*1 





a 

-d 






-d 
a 







b 



(19) 



*2*3 + *0*1*2*3 
*2*3 — *0*1*2*3 



(20) 
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and 



(^0)2,3 



--1 



/ -a' 

-b' 

d 

V d! 







V 




-d>\ 




a' J 



2,3 



/ -v 




-d 



\ 



d \ 

-a' d! 

-d! a! 

V ) 



where 



a = 



c = 



-hh - *o*3 + i*o*2 + 2*1*3 
-*i*2 + *o*3 + itote - 2*1*3 



b' = -ht 2 - t t3 - it t2 - itih , 
d! = -tit 2 + t t3 - it t2 + itih . 



(21) 



(22) 



We regard these matrices as two by two block matrices with 2-dimensional blocks. The indices 
(«2, J2) numerate the blocks while the indices (^3,^3) numerate matrix elements inside the blocks. 
From (7) it follows that 



ab = a'b' , 



cd 



dd', 



2 1 i2 2 

a + b — c 



d 2 = 0, 



a' 2 + b> 2 



d 2 -d' 2 



0, (23) 



therefore each of the above matrices can be viewed as an i?-matrix of the eight-vertex free- 
fermion model, which satisfy the condition 



^(FFM) 



/ uj\ ujj ^ 

W3 W5 

ujq uji 

\ OJg U>2 ) 



(24) 



3.2 Elliptic parameterisation of the weights 

Here we will show that the weights (10) can be parametrised in terms of Jacobi elliptic functions. 
First, note that the form of the matrices (19) and (21) and the relations (23) suggest to use 
Baxter's parameterisation of the eight-vertex model [33] (the latter contains the symmetric 
free-fermion model as a particular case). To do this define the elliptic modulus 



k 



cd 
ab 



ddf_ 
a'b' ' 



(25) 



similarly to that in the symmetric eight-vertex model with the weights a,b,c,d. (cf. Eq.(5.7) 
of [33]). With an account of (20) and (22) the RHS of Eq.(25) is a function of the dihedral 
angles 61,62,03. After elementary simplifications, one obtains 



k 



1 — sin < 



sin ( 



4*0*1*2*3 



l + sin^' (1-*2*?)(1 -*|*§) 

Then using the relations (see §132 of [34]) 

a 3 



oi7 . ai . «2 
sin an = IK sin — sin — sin ■ 

2 2 2' 



sin ctj 



1 V sin «o sin at\ sin a 2 sin 03 
sin#i 



2K sin — cos — cos — , 
2 2 2 ' 



where (i, j, k) is any permutation of (1, 2, 3) and K is defined in (3), one obtains 

2 sinai . . ■ a ■ 

sin <p = K — sin 02 sin 03 = sin 62 sin a 3 

sin 6\ 



(26) 



(27) 



(28) 



which shows that (f> is the angle between the line 1 and the plane containing the lines 2 and 3, 
as illustrated in Fig.3. Note that value k = corresponds to the right angle (f> = tt/2. 
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Figure 3: This figure clarifies the geometric interpretation for the angle </>, given by (28), as an angle 
between the line 1 and the plane containing the lines 2 and 3. 



Further, let sn x = sn(x,k), en x = cn(x,k), dn x = dn(x, k) and cd x = cd(x,k) denote 
Jacobi elliptic functions of the argument x and modulus k (we follow the notations of [35]). Now 
let us parametrise the angles # 2 and #3 in terms of new parameters wi and u> 2 , 

e -i8 2 = (l+sn(2 W i))(l-fc S n(2^i)) ^ = (1 + snQ; 2 ))(l - fc sn(2w 2 )) 

cn(2u>i)dn(2^i) ' cn(2w 2 ) dn(2w 2 ) ' 1 ' 

Given 62,63 each of these equations have two solutions in the periodicity rectangle (we assume 
< k < 1) 

< w 1 ,w 2 < 2K, < Im Wl ,Imw 2 < iK' /2 . (30) 

where K = K(/c) and K 7 = K(\/l — k 2 ) are the complete elliptic integral of the first kind. Now 
substitute (29) back into (26) and solve the resulting equation for 6\ in terms k,w\,W2- There 
are four solutions of the form 

= ^4 ■ (31) 
cd(2i«i) v ; 

corresponding to four possible choices of w\,W2 in (29). Only one of them 3 lead to the original 
value of 6\. This allows one to uniquely fix the required solution of (29). 

From now on we regard k,wi,W2 as new independent variables instead of 61,62,63. The 
latter are now defined by the formulae (29) and (31). The weights in Eq.(10) are expressed 
through the square roots of tangents of halves of spherical excesses ao, «i, 012, 013. Surprisingly 
enough, after lengthy calculations we found that these tangents have rather simple expressions 
in the new variables, 



1 . 1 

2 



to = (-if(w-)/f(K-w+)y, h = (»/(«;_) /(K -«/+)) 

t 2 = (-if(K-w„)f(w + )y, t 3 = (i/(K -«;_)//(«;+)) 



(32) 



Calculating 6\ from (31) for all four solutions of (29) one obtains the following values 8\ = ±d[ true \ ±9[ 
where o[ true ^ is the original value of 81, while 6[ is a completely different value, corresponding to a different 
spherical triangle, (#i,#2,#3) which, nevertheless, leads to the same modulus k in equation (26). 
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where 

w± = wi ± w 2 , k' = y/l — k 2 
The function f(w) is defined as 

v k sure 

/H 



(cnro + dnw)(& cnw + dnw) 
It follows that the products appearing in (20) and (22) simplify to 

t h = f(w-), t 2 t 3 = f{K-w. 



(33) 



(34) 



( f(w-)f(K-w-) 
03 " \f(w + )f(K- W+ ) 



f(w + ) f(K - w+)t Q t 3 . 



(35) 



Then using addition theorems for elliptic functions, one obtains 

a = p- cdw_ , b = p_ snw_ , c = p- , d = p- k cdw_ snw_ , 
a' = p + cd w + , b' = p + sn w + , d = p+ , = p+ fc cd w + sn w + 

where 

4(1 — snro_) dn w_ 



(36) 



P- = 1 + *0*1 + *2*3 - *0*1*2*3 = 

P+ = t t3 ~ hh + i{tot 2 - ht 3 ) = 



cnw_ (cnw_ dnw_ + (1 — snw_)(l + k snw_)) 
'cdw_ snio_ 



^cdw+ snw +< 

Next, we specify the parameter £ in (17), which so far remained at our disposal, 



£ = -y- = h(wi,w 2 ), h(x,y) 



' cd{x — y) sn(x — y) 
cd(x + y) sn(x + y) 



(37) 



(38) 



Finally, introduce the .R-matrix of the symmetric eight-vertex model [33] specialized to the 
free-fermion case, 



H = 



/ cdw k cdw snw \ 

snw 1 

1 snw 

\ cdw snw cdw / 



(39) 



Using this notation, the matrices (19), (21) can be written in a uniform way 

(L°) 23 = p_R 23 (w_), (Lj)23 = -p + r 1 4 2) M 23 (w+), 

(L?) 23 = i P+ ia^R 23 (w + )a^\ (L\) 23 = 



(2) 



w- a, 



(2) 

y i 



(40) 



(2) (2) (2) 

where the Pauli matrices a x , o~ y , o~\ act in the space 2. Taking into account (37) and (38) is 
easy to see that all expressions for the matrix elements in (40) are meromorphic double periodic 
functions of the variables w\ and w 2 . 

To summarize, we have shown that all matrix elements the operator 

Ci 23 { Wl ,w 2 | k) = Li23(0i, 2) ft) = (D(£) <g> F ® F) Ri 23 (0i, 3 ) (D(£) ®F0 F)" 1 , (41) 

are meromorphic functions of w\ and w 2 , provided the two sets of variables, entering different 
sides of this equation, are related by (26), (29) and (31) and the parameter £ is defined by (38). 
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We would like to stress the above elliptic parametrisation was obtained for a generic case 
without any special requirements for the values of 01,02,03- By construction, this parametri- 
sation breaks the symmetry of the weights (the directions 1 is distinguished). Moreover, the 
modulus k will, a priori, be different for different vertices of the tetrahedron. Therefore, in 
general, one cannot parametrise all weights in the tetrahedron equation by elliptic functions 
of the same modulus. Nonetheless, there exists an important four-parameter reduction of this 
equation, where such parametrisation is possible. This is the "prismatic limit" considered in the 
next section. 

4 Reductions of the tetrahedron equation 

In this section we consider certain reductions of the tetrahedron equation for the solution (10) 
and its connection to an example of the tetrahedral Zamolodchikov algebra constructed in [19] . 

4.1 The "static limit" 

The term "static limit" stems from the original Zamolodchikov's work [1] where it was related 
to the case of slowly moving or "non-relativistic" straight strings in 2 + 1 dimensions. On the 
level of parameters this corresponds to a configuration where the sum of dihedral angles is equal 
to 7r for every vertex of the tetrahedron. Namely, for Eq.(15) this means 

01 + 02 + 03 = 7T, 01 + 04 + 05 = IT, 

(42) 

04 + 06 - 02 = 0, 03 + 06 - 05 = 0. 

The linear angles between edges at the vertices all become equal to or ir and the tetrahedron 
degenerates into four planes intersecting along the same line. Their relative orientation is fixed by 
three angles only, so the tetrahedron equation in this case contains three independent parameters 
(instead of five in the general case). 

We will use a special notation for the operator R(0i, 02, 03) specialized to this case 

S(0i,0 2 ,0 3 ) = R(0i,0 2 ,0 3 ), # 1 + 2 + # 3 = 7r. (43) 

Its matrix elements are determined by (10) where one sets 



t = 0, U = Ti = ytan^, i = 1,2,3. (44) 
There are only sixteen non-vanishing elements, 



o0,0,0 
"-'0,0,0 


_ c-0,1,1 _ 
- ^0,1,1 - 


^1,0,1 - 


c+1,0 _ 

^1,1,0 - 


1 


q0,1,0 
0,0,1 


c0,0,l 
— °0,1,0 — 


ol,0,0 _ 


J i,o,o — 


T 2 T 3 


c+1,1 
D 0,l,0 


_ c+0,0 _ 
- °0,Q,1 ~ 


q0,0,1 
^1,0,0 - 


o0,l,0 _ 
-^1,1,1 - 




c0,0,l 

1,1,1 


r.1,1,1 


rrl,0,0 

^Ojl.O — 


q0,1,0 

'-'i.o.o — 


TiT 2 



(45) 



This solution of the tetrahedron equation (with a slightly different parameterisation, see Sect. 4. 5 
below) was first obtained in [19]. Two years later [17] it was understood as the vertex form of 
the original Zamolodchikov's solution in the static limit [1]. 
The operator S satisfies the relation 

S 2 = 1. (46) 
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It has a block-diagonal form with two 4x4 blocks, a trivial one, which coincides with the identity 
matrix, and a non-trivial one which has two eigenvalues +1 and two eigenvalues —1. 
The operator S possesses left and right "bare vacuum" eigenvectors 

(vo v vq) 1 S = (vq vq vqY , S (v v v ) = i>o «o v , v = . (47) 
Note that this property does not hold for the general case (10) with to / 0. 
4.2 The planar limit 

Geometrically the "planar limit" corresponds to the case when all four vertices of the tetrahedron 
lie in one plane, i.e. when the tetrahedron becomes "squashed" into a plane. An edge-spin 
solution of (13) corresponding to this case was first obtained by Hietarinta [20]. Subsequently it 
was understood as the planar limit [36] of the Zamolodchikov model. The details of this reduction 
are rather complicated and not immediately related to the main topic of this paper. Therefore, 
we refer interested readers to the original publication [36] where this limit is thoroughly studied. 



4.3 The prismatic limit 

Here we consider yet another limiting case of the equation (13), when one of the tetrahedron 
vertices goes to infinity (we choose it to be the vertex corresponding to Ri23). Then the sum of 
dihedral angles at this vertex will satisfy an additional constraint, 

01 + 2 + 03 = it, (48) 

whereas the whole tetrahedron turns into an infinite prism. The number of independent angles 
reduces from five to four. The edges 1, 2 and 3 become parallel and therefore, have the same 
angle to the plane, containing the edges 4, 5 and 6, forming the base of the prism. Remembering 
the geometric definition (26) of the elliptic modulus k in the previous section, we conclude that 
all weights corresponding to the vertices (1,4,5) , (2,4,6) and (3,5,6) can be parametrised by 
elliptic functions of the same modulus. 

Eqs. (26), (29) and (31) define a generic change of variables from three dihedral angles 
{01,02,0s) to new parameters (k, w±, it^)- Here we want to apply this substitution for three 
different sets of dihedral angles corresponding to the three vertices at the base of the prism 

(01,04,05) -> (A:, 1*1,1*2); (?r- 02,04,0e) -> (A;, 1*1,1*3); {0 3 ,ir - 5 ,0 e ) ^ (k,u 2 ,u 3 ). (49) 

As explained above from the geometric considerations the elliptic modulus k will automatically 
be the same for all three sets. For instance, let (bi, 62, H) be the sides of the spherical triangle 
with angles {0i,04, #5), then we set 

k = i — %. sin 6 = sin 64 sin 63 . (50) 

1 + sin <p 

Next define the variables u±,U2, 1*3 by the relations 

e _ m = (l + S n(2m))(l-fcsn(2m)) 

cn(2i*i) dn(2iti) 
^ = (l + sn(2n 2 ))(l-fcsn(2n 2 )) 

cn(2tt 2 ) dn(2it 2 ) 
= (l + sn(2H 3 ))(l-fcsn(2n 3 )) 

cn(2tt3) dn(2it3) 

and 

m = cd(2u 2 ) eie2 = cd(2ui) ei03 = cd(2n 3 ) 
cd(2ui)' cd(2u 3 )' cd(2n 2 ) 



(51) 



(52) 



11 



which are obtained by a simple specialisation of (29) and (31) for the three substitutions (49). 
Note that the above formulae (51) and (52) give a parameterisation of six angles of a triangular 
prism in terms of elliptic functions. In particular, it is easy to check that the three expressions 
in (52) are consistent with (48). 

Let us now rewrite the tetrahedron equation (13) for this case in the new variables. We will 
do this in several steps. First, guided by the formula (41) define three operators 

£i45(«i,«2) = (D(£i)<g>F<g>F) Ri45(0i, #4 A) (D(6)®F®F)-\ 

£246(«1,«3) = (D(£ 2 ) <g> F <g> F) R 2 46(*-- 02,04,0 6 ) (D(&) <8> F <g> F) _1 , (53) 

£356(u 2 ,u 3 ) = (D(£ 3 )®F®F) R 145 (03,^- 05,0 6 ) (D(e 3 ) <E> F (8) F)- 1 , 

Here F and D(£) are defined in (18) and the parameters £1, £2, £3 have been chosen in agreement 
with (38) in each of the three cases, 

€i = h(ui,U2), & = h(u!,u 3 ), & = h(u 2 ,u 3 ), (54) 

where h(x,y) is defined in (38). Therefore, according to the result of Sect. 3, all operators (53) 
are meromorphic double periodic functions of u±,u 2 ,u 3 , which, of course, implicitly depend on 
the elliptic modulus k. 

Next, we need to express in the new variables the weights, corresponding to the infinitely 
distant vertex (1,2,3). It is convenient to define a new operator, 

5i 23 («i,«2, « 3 ) = (D(6) D(6) ® D(&)) S 123 (0i, 2 , 3 ) (D(6) D(&) Dfe))- 1 , (55) 

which differs from (43) by a diagonal similarity transformation. 
Using (52) it is not difficult to show that 

tan ^ = -ig(m, u 2 ), tan ^ = — — r, tan ^ = -ig(u 2 , u 3 ), (56) 

z l g{ui,u 3 ) I 

where 

, 2 sd(x - y) sd(x + y) 

9{x,y) = k — — - — — , (57) 

cn(x — y) cn(x + y) 

Combining this with (45), one obtains the following expressions for the non- vanishing matrix 
elements of S(ui,u 2 , u 3 ) 

cOOO _ cllO _ clOl _ cOll _ i (zo\ 

50%° = U(u u u 2 ,u 3 ), = U(u 3 ,u 2 , Ul ), S^ = V(ui,u 2 ,u 3 ), 

5 ioo = U(ui,-u 2 ,-u 3 ), S$fi = U{u 3 ,-u 2 ,-ui), = V(-u u u 2 ,-u 3 ), 

S™{ = U( Ul ,-u 2 ,u 3 ), cSin = U{u 3 ,-u 2 , Ul ), = V( Ul ,u 2 ,-u 3 ), 

5 doi = U(u!,u 2 ,-u 3 ), cS^q 1 = U(u 3 ,u 2 ,-ui), = V(-U!,u 2 ,u 3 ), 

where 

, . sd(x + y) cn(x - z) /2 sd(y - x) sd(y + z) 

U(x,y,z) = — — - — ■ — -, V(x,y,z) = —k — — ■ — r -, 59 

cn(x — y) sd(x + z) cn(y + x) cn(y — z) 

It is easy to see that all above matrix elements are meromorphic functions of u\,u 2 ,u 3 , implicitly 
depending on the elliptic modulus k. Note that even though the parameters £1,^2, £3, entering 
(55), were fixed by the analyticity requirements for the other three operators, defined in (53), 
the matrix elements of S(u±,u 2 ,u 3 ) have automatically became meromorphic functions on the 
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elliptic curve 4 . Note also that the similarity transformation in (55) introduces two additional 
parameters into (58) and breaks the symmetry relations between different matrix elements exhib- 
ited in (45). Remind that the original static limit weights (45) depend on only two independent 
angles. 

To complete our analysis of the prismatic limit apply the similarity transformation with the 
matrix 

G = D(£i)®D(£ 2 )®D(£ 3 )®F®F®F (60) 

to both sides of the tetrahedron equation (13). Then using the definitions (53) and (55) one 
obtains 

5l23 (u\ ,U 2 ,U 3 ) £145 (ui , u 2 ) C 2m (ui , lt 3 ) £ 356 (U2 , u 3 ) 

(61) 

= £356(U2,U 3 ) £ 2 46(U1,U 3 ) £u 5 (ui,U2)S 12 3{u 1 ,U 2 ,U 3 ) . 

As required, this equation contains exactly four independent parameters, namely, 1*1,1*2,1x3, 
which are indicated explicitly, and the elliptic modulus k which is implicitly assumed. 



4.4 Elliptic parameterisation of the static limit 

Consider the tetrahedron equation (13) in the static limit, i.e., when the angles satisfy the 
constraints (42). In the previous subsection we saw that upon the change of variables (52) and 
the similarity transformation (55) the vertex weights £123(1x1, 1x2, 1*3), corresponding to the vertex 
(1,2,3), become meromorphic functions on the elliptic curve. Remarkably, it is possible to do 
this for the other three vertices as well, using the elliptic functions of the same modulus. First, 
we need to parameterise the angles, satisfying (42), in term of the elliptic functions. To avoid 
confusions with the notations of the previous subsection we denote these angles 9i, 6 2 , . . . , 9q. 
Let us parameterise 9\,9 2 ,9 3 by the same formulae as (52), 

is-, = cd{2u 2 ) j- d2 = cd(2m) ih = cd(2ix 3 ) , . 

cd(2ixi)' cd(2n 3 )' cd(2ix 2 ) 1 ' 



and assume that 

g ie 4 = cd(2m) # B = cd(2n 4 ) = cd(2?x 4 ) 

cd(2i* 4 )' cd(2n 2 )' cd(2u 3 ) 

This implies the expressions (56) with 9±,9 2 ,9 3 replaced by 9\,9 2 ,9 3 and 



(63) 



tan^ = — — -, tan^ = -15(1x2,1x4), tan ^ = -ig(u 3 , 1*4), (64) 

2 g{ui,U4) 2 2 

where g(u, w) is defined in (57). Note, that these formulae contains five independent parameters: 
1x1, ix 2 ,ix 3 , 1x4 and k, whereas there are only three independent angles, satisfying (42). So the 
above parameterisation contains two spurious parameters. 
Next, assume that £1,^2, £3 are given by (54) and 

£4 = /i(ui,U 4 ), £5 = /*(li 2 ,lt4), £ 6 = /t(lX 3 ,l*4), (65) 

where h(u,w) is defined in (38). Applying now the similarity transformation with the matrix 

G = D(&) ® D(6) ® D(&) ® D(£ 4 ) ® D(£ 5 ) D(£ 6 ) (66) 

to both sides of the tetrahedron equation (13) and taking into account the above parameterisa- 
tion, one obtains 

Si2 3 (ui,u 2 ,u 3 ) 5 245 (1x1 , 1x2 , 1x4) 5 246 (1x1 , ix 3 , 1x4) 5 356 (u 2 , ix 3 , 1x4) 

(67) 

= <S 356 {u 2 , IX3 , IX4) 5 24 6 (iXi , IX3 , IX4) 5 245 (iXi , 1*2 , IX4) 5i 23 (U\ , 1*2 , 1*3 ) 



4 Eq.(61) below implies that ratios of matrix elements of 5i23(wi, M2, W3) are meromorphic functions on the 
elliptic curve. However, some of these matrix elements in (58) are equal to one, therefore the remaining elements 
should be meromorphic. 
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where the matrix elements of Sijk(u,v,w) are defined by (58). It appears, that this equation 
naturally complements the equation (61). As noted in Sect. 4.1 there are only three independent 
angles in the tetrahedron in the static limit. However, Eq.(67) contains five independent param- 
eters (four u's and the modulus k). Two additional parameters are not spurious in this case, 
as they cannot be removed by a re-parameterisation. They were introduced via the similarity 
transformation with the matrix (66). 

4.5 Tetrahedral Zamolodchikov algebra 

Here we show that some additional reduction of the tetrahedron equation (61) leads precisely 
to the example of the tetrahedral Zamolodchikov algebra constructed in [19]. It is convenient to 
introduce two 4x4 matrices, acting in the product of two spaces C 2 (g) C 2 as 

n a (u,w)= £ Civ!; (<<■>>■ O a = o,i, (68) 

»2,»3,j2,j3 

where (Ef)^ = Si a 5jb denotes the 2x2 matrix unit. According to (40), (41) these new matrices 
are simply related to the i?-matrix of the free-fermion model (39), 

n°(u,w) = p-R(u-w), K 1 {u,w) = -p + C 1 (vz®l)^(u + w). (69) 

Now calculate matrix elements of tetrahedron equation (61) sandwiched between the fixed vec- 
tors 

(L\ = (vl ® vl <g> v£)i2 3 , \R) = (v a ®v b ®v c ) 123 (70) 

in the first three spaces, where the basis vectors vo, v\ are defined (2) and the indices a, b, c take 
the values 0, 1 (in the notations of (15) this corresponds to the case when i\ = i 2 = i 3 = and 
ji = a, 22 = b, j 3 = c). The resulting equation considerably simplifies, since the matrix S in the 
LHS of (61) drops out due to (47). Writing this equation in an abbreviated form, one gets 

(£g) 4 5 (£o)46 (^8)56 = ^2 ^def (£o)56 (£o)46 (£o)45 (71) 

d,e,f 

where, for example 

(£g) 45 = (uo ® 1 ® 1) £145 (v a <8> 1 <8> 1) (72) 

and similarly for the other C's. For fixed values of a, 6, c, both sides of (71) still remain operators 
acting in the spaces 4,5,6. At this point it is convenient to relabel these spaces as 1,2,3 (the 
former spaces 1,2,3 are no longer required, so there should be no confusions) and use the 
definition (68). In this way one obtains 

1Z\ 2 («i , u 2 ) Tl\ z (tii , u 3 ) K% 3 (u 2 ,u 3 ) 

= 2 S2 b e c f (u 1 ,U2,u 3 ) K f 23 {u 2 , u 3 )n e 13 ( Ul , U3 ) Ki 2 ( Ul ,u 2 ) 

d,e,f 

This relation is known as the definition the tetrahedral Zamolodchikov algebra. It was intro- 
duced in [19]. The coefficients S^j(ui,u 2 ,u 3 ) can be considered as the structure constants 
of the algebra. They satisfy the tetrahedron equations (67), which play the role of the asso- 
ciativity condition for the defining relations (73). This relationship is analogues to that for 
the Zamolodchikov-Faddeev algebra [37] where the Yang-Baxter equation plays the role of the 
associativity condition. 

In [19] Korepanov constructed a representation of the the algebra (73). He essentially pos- 
tulated the expressions (68), without any connection to the 3D Zamolodchikov model and then 
calculated the coefficients S^f(ui,u 2 ,u 3 ), solving the linear equations (73). He used the nor- 
malization, where £1 = £2 = £3 = 1, so his expression for S123 (given before Eq.(2.32) in [19]) 
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(73) 



contains square roots and, upon a change of variables, exactly coincides with our expression 
for S(#i, #2, #3), given by (45). More precisely, Korepanov used the variables Ai,A2,A3 (which 
exacly correspond to our variables iti,it2,it3, respectively) and the variables (pi,<f2,f3, defined 
by Eq.(2.32) of [19], 

1 — cd(2m) 

tanhy,= * , * = 1,2,3. 74 

1 + cd{2ui) 

Using addition theorems for elliptic and trigonometric functions it is not difficult to show that 

g(v,j,u k ) =tanh(tpj -tpk), j, k = 1,2, 3. (75) 

where g(x,y) is given by (57). Then it follows from (56) that 

61 = 2i(<p 2 - <pi), 62 = vr + 2i(<pi - </? 3 ), 63 = 2i(<p 3 - Lp 2 ) . (76) 

Subsequently, Shiroishi & Wadati [22] reproduced the calculations of [19]. They worked in 
essentially the same normalization and notations as this work and their expressions for the matrix 
elements of £123 (given by their Eq.(5.4)) exactly coincide with our Eq.(58) with interchanged 
u\ and U2- In addition, Eq.(73) was re-checked once again in [23] in the trigonometric limit 
(fc = 0), see §12. A. 1 in [23]. 



5 Conclusion 

In this paper we have considered a special limit of the tetrahedron equation in the 3D Zamolod- 
chikov model, where one of its vertices goes to infinity and the tetrahedron turns into an infinite 
prism. We explicitly showed that all vertex weights in the tetrahedron equation (61) in this case 
are meromorphic functions depending on three points on an elliptic curve. The parameterisation 
of the weights are given by (40) and (58). As a byproduct we obtained a parameterisation of 
angles of a triangular prism in terms of elliptic functions (a problem, which is, certainly, too hard 
for school geometry, but could have taken a comfortable place among exercises in the classic 
textbook by Whittaker & Watson [35]). The corresponding formulae are given in (51) and (52). 

Next, we have shown that a further reduction of this, already special, case of the tetrahedron 
equation leads precisely to the "tetrahedral Zamolodchikov algebra" , originally invented in [19] 
without any connection to the 3D Zamolodchikov model. Interestingly, this algebra (more 
precisely, its trigonometric variant, corresponding to k = 0) has been used by Shiroishi & 
Wadati [21] for a two-layer construction of the Shastry's solution [24] of the Yang-Baxter equation 
connected to the Hubbard model [38]. Recently, the interest to this solution was renewed due 
its remarkable appearance [25-28] in the problem of the AdS / CFT correspondence for M = 4 
SUSY Yang-Mills theory in four dimensions. 

Actually, one of motivations for this work was to unravel an algebraic origin of the Shastry's 
i?-matrix, which so far has not been included into the standard quantum group scheme [39]. 
This .R-matrix has long been suspected to have a hidden 3D structure. Indeed, it has an evident 
two-layer structure implied by the construction of [21], mentioned above. Moreover, it does not 
possess the "difference property", which is another indication of a 3D origin. To enhance the 
last argument it is worth mentioning that the chiral Potts model [40,41], which is the most 
notable example among solved models without the difference property, is just a two-layer case 
of the multi-state generalisation of the 3D Zamolodchikov model [7,42,43]. 

Our results suggest that the Shastry's i?-matrix is closely related to the 3D Zamolodchikov 
model, which was originally formulated as an exact relativistic factorized scattering theory of 
"straight strings" in 2+1 dimensional space-time. This 3D model has an exremely rich algebrac 
structure, previously studied in [5-9,13-17,32]. It would be interesting to see how this connection 
will help to understand real origins of the Shastry's i?-matrix. We postpone these considerations 
to a separate publication [44]. 
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